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Favorable growing conditions for fine 
fuels in Hawaii sustain an environment 
prone to conflagrations. In 1973, the State 
of Hawaii began use of Fuel Model A of the 
1972 National Fire-Danger Rating System 
to diagnose fire potential. Soon after, a 
study was begun of the physical character- 
istics of three major fine fuels: Andropogon 
virginicus (broomsedge), Melinis minuti- 
flora (molasses grass), and Pennisetum 
setaceum (fountaingrass). Surface area-to- 
volume ratios were estimated with the 
reciprocal harmonic mean. An approxima- 
tion of the complete elliptic integral of the 
second kind was used as an additional 
refinement to surface area-to-volume ratio 
estimates of stalks with elliptical cross- 
sections. Bed depths and fuel loadings were 
estimated, also. Loading estimates were 
apportioned into live and dead categories 
on the basis of observed proportions-by- 
weight of each class. Surface area-to- 
volume, loading, and fuel depth estimates 
were generally higher than the corre- 
sponding parameters of standard grass 
models. 

Retrieval Terms: Hawaii fine fuels, surface 
area-to-volume ratio, fire-danger rating, 
harmonic mean 


F avorable growing conditions for 
fuels aggravate the problem of 
wildfire in Hawaii. A forest ravaged 
by fire soon gives way to grasses that 
threaten to hasten the spread of the 
next wildfire. Examples of this cycle 
are found throughout the State — at 
Makua and Helemano on the island 
of Oahu, Puna on the island of 
Hawaii, Puu Ka Pele on Kauai, 
Kahikinui on Maui, and Kaunakakai 
on Molokai. 


To diagnose the day-to-day fire 
potential on conservation land, the 
Hawaii Division of Forestry began 
using Fuel Model A of the National 
Fire-Danger Rating System, in 1973, 
on a limited scale. ^ The Division sub- 
sequently phased into the AFFIRMS 
computerized system for processing 
data to estimate fire danger. 2 In the 
meantime, we started a study of the 
physical characteristics of three major 
fine fuels in Hawaii: broomsedge 



Figure 1— Broomsedge was sampled on the fringe of the Helemano (Poamoho) Forest 
Reserve, Oahu; molasses grass In the Kuaokala Forest Reserve near Ka’ena, Oahu; and 
fountaingrass in the Keamuku area of South Kohala, island of Hawaii. 



(Andropogon virginicus), molasses 
grass (Me Unis m inn! (flora), and 
fountaingrass (Pennisemm setaceum). 

This note reports estimates of 
loading, surface area-to-volume ratio, 
and bed depth of the three grasses 
studied. The results suggest that 
standard fuel model parameters 
generally underestimate the character- 
istics of Hawaii fine fuels. 


STUDY SITES 

We sampled broornsedge on the 
fringe of the Helemano (Poamoho) 
Forest Reserve, molasses grass in the 
Kuaokala Forest Reserve near ICa’ena, 
and fountaingrass in the Keamuku 
area of South Kohala. These fire- 
prone areas offered pure stands of the 
respective fuel types (Jig. ]). 

The broornsedge site at Helemano 
is on a modest slope (locally fiat) west 
of the Koolau Range, at 427 m 
elevation. Median annual rainfall at 
the Helemano Forest Reserve for 40 
years before 1958 was 207 cm.-' In 
1975, annual rainfall was 149 cm.”^ The 
site of molasses grass sampling is 
leeward of the Waianae Range ridge- 
line, and considerably drier. At the 
Kuaokala Forest Reserve, the rain 
gage nearest the sampled plots 
recorded a median annual rainfall of 
81 cm in a 28-year period. In 1975, it 
recorded 75 cm. M.olasses grass was 
sampled at 442-m elevation. Fountain- 
grass was sampled in the area where 
the State experiences its largest wdld- 
fires (fig, 2). In winter 1978, a fire 
swept through 12,000 acres (4856 ha) 
of Keamuku in Ifi days. Prolonged 
dry spells and strong, persistent winds 
can result in extremely dangerous fire 
conditions in the fountaingrass stands 
that blanket the vast lava fields, 
sloping gently upward from the Kona 
coast. Annual rainfall in the Keamuku 
area was 67 cm in 1975, compared to 
the median for 56 years of 71 cm. Site 
elevation is about 762 m. 


re 2— Three major fine fuels studied 
(A) broornsedge, sampled at Hele- 
D, Oahu; (B) molasses grass, sampled 
(a’ena, Oahu; and (C) fountaingrass, 
)led at Keamuku, island of Hawaii. 





EXPERIMENTAL DESIGN 

The sampling design was largely 
that of Brown, 5 and BurganT We 
sampled in the fuel bed at 10 points 
spaced 50 feet (15.24 m) apart, on a 
randomly placed linear transect. Four 
subplots were established at each 
point {fig. 3). Each subplot measured 
40 X 80 cm, and was oriented parallel 
to the transect. By using this design, a 
total of 40 subplots per transect w^ere 
sampled. 

Within each subplot, we measured 
depths of leaves and stalks from each 
of three points at 20-cm intervals 
along the midlength of the subplot. 
Depth corresponded to the highest 
fuel particle intersecting an imaginary 
4-inch (lO-cm) diameter cylinder 
centered on the depth sampling point. 
Anderson'7 suggests that this depth be 
multiplied by 0.7 to obtain the best 
estimate of depth that influences fire 
behavior. We singled out the subplot 
with the greatest volume of plant 
material — the base subplot — and esti- 
mated the volume of material in the 
other subplots, relative to the base 
subplot. Because we depended on 
observer skill, we decided to calibrate 
observer judgment in field training. 
The volume estimates in this report, 
therefore, are the average of two 
independent observations, which 


generally agreed within 10 percent of 
each other. 

To estimate fuel loading, we needed 
to remove the material from the 
subplot. We did this only on the base 
subplot. A portable frame about 5 feet 
(1.52 m) tall, constructed to the length 
and width of the subplot, was used to 
define the subplot volume. Stalks and 
leaves protruding from the frame were 
trimmed back to the face. After the 
volume was so defined, the materia! 
was clipped as near the ground as 
possible and sealed in plastic bags. 
The samples were air-dried and 
weighed. 

Additionally, observers estimated 
the proportion of live (green) material 
in the base plot. Again, some field 
training was required to ensure reli- 
ability of observer judgment. It was 
generally more difficult to estimate 
this proportion than to estimate rela- 
tive volume with respect to the base 
plot. We obtained an objective 
proportion-by-weight by separately 
w'eighing green and dead material in 
broomsedge and fountaingrass. 

To estimate the surface area-to- 
volume ratios (cj), the simplest as- 
sumptions were made about stalk and 
leaf shapes. We assume that stalks are 
right cylindrical, if we denote surface 
area by S and volume by V, and 
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Figure 3— Sampling was done in the fuel bed at 10 points spaced 50 ft (15.24 m) apart, 
on a randomly placed linear transect. Four subplots were established at each point, 


assume that the end areas of the 
cylinder are negligible, then 


S 2m I 

r d 
in which 


Eq. I 


r radius of cylinder 
/ = length of cylinder 
d = diameter of cylinder 

For a simple leaf, let t be leaf thick- 
ness. Ignoring the area on the edge, 



2 


Eq. 2 


Brown^ discusses calculating surface 
area-to-voliime ratios by using photo- 
micrography. Estimates are obtained 
from Equations I and 2, with the 
arithmetic mean of stalk diameters 
substituted for d and the arithmetic 
mean of leaf thicknesses substituted 
for t. We choose to use the harmonic 
mean instead, for these reasons: 

Assume that the arithmetic mean is 
the ‘'•'besr estimator of some quantity.^ 
Suppose we have n independent leaf 
thickness measurements (x,, X 2 , . . . x„). 
When we compute the arithmetic 
mean of thickness 


X is the best estimate of leaf thickness 
we obtain from our sample. 

We want the best estimate of 
surface area-to-volume ratio, however. 
The assumptions of the preceding 
paragraph lead to the conclusion that 
surface area-to-volume ratio is a 
function of the inverse of thickness. 
Our same sample of n leaf thicknesses 
can then be considered a sample of n 
independent surface area-to-volume 
ratios (2/Xj, 2 /X 2 , ... 2/xJ. It 
follows that the best estimate of 
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surface area-to-volume ratio is the 
arithmetic mean of the {2/X|]: 


^leaf 



n ^ X: 
i = ! 


The same logic applies to the best 
estimate of surface area-to-volume 
ratio for stalks, with the more general 
result 


c 


Z - 

n ^ X: 


Eq. 3 


they were less elliptic than the foun- 
taingrass stalks. 

By assuming a right elliptic cylin- 
drical fuel, the surface area-to-volume 
calculation is encumbered by an 
elliptic integral, namely, a complete 
elliptic integral of the second kind, 
which determines the perimeter of the 
fuel cross-section normal to the length. 
We use an approximation of the 
elliptic integral, and obtain 





vr(f)(|)/ 


we can express Equation 4 as 


o 


4 { y [ c 2 + c 2 ( 1 - 62)]}'/2 



c[cMI-e2)] 


Vi 


The approximation to the elliptic 
integral, 



in which c is 2 for leaves and 4 for 
stalks. The summation is the recip- 
rocal of the harmonic mean, H, of 
thickness (diameter) measurements. 
Denote the arithmetic mean of thick- 
ness (diameter) by A; it can be shown 
by the Cauchy-Schwartz Inequality 
for Sums^o that for finite Xj > 0, I / H 
> 1/ A; equality obtains when the {xj} 
are all equal. By using the arithmetic 
mean of thickness measurements in 
Equations 1 and 2, therefore, we 
underestimate o. We should, instead, 
use the reciprocal harmonic mean, as 
in Equation 3. 

We randomly measured leaf thick- 
nesses and stalk diameters on the air- 
dried samples collected from the base 
subplots. Samples were scattered 
evenly on a table. By using the table 
edges as coordinate axes, and a table 
of random numbers, we generated 
bivariate random coordinates, at 
which we located a leaf or stalk 
sample. A vernier caliper with 0.001- 
inch graduations measured thickness 
and diameter. Twenty measurements 
were obtained for each subplot on 
even-numbered points in the transect. 
The surface area-to-volume ratio was 
estimated on the basis of a sample size 
of 100. 

DISCUSSION 

The stalk cross-sections, particu- 
larly in fountaingrass, were sometimes 
more elliptic than circular. We at- 
tempted, therefore, to measure the 
major and minor axes of the ellipse on 
fountaingrass stalks. We did the same 
on molasses grass stalks, although 


4 [ y ( c 2 + d 2)]'/2 
= Eq.4 

in which 

c = major axis 
d = minor axis 

When c = d. Equation 4 becomes 
identical to Equation 1. By definition 
of eccentricity, e, of an ellipse. 


- 7r[y (c^ + d^)]/^ 


which first appears in Equation 4, is 
good, even for high eccentricity. With 
e = 0.7071, for example, the error of 
the approximation is 0.0072 (on the 
basis of 5 significant digits). 

The statistics for stalk and leaf 
surface area-to-volume ratios were 
computed from the inverse of thick- 
ness measurements. They were gener- 
ated by BMDP Program 2D,‘2 which 
computes descriptive statistics for a 


Table \ — Minimum and maximum, mean, median, variance, and standard error of the mean for stalk and 
leaf depths of three major grass fuels in Hawaii 


Species 


Min 

Max 

Mean 

Median 

Variance 

Standard 
error of 
the mean 













0 (0) 

41 (104) 

26.86 (68.22) 

31.0 (78.7) 

107.33 

(692.45) 

0.95 (2.41) 

3 (8) 

28 (71) 

15.74 (39.98) 

16.0 (40.6) 

19.03 

(122.77) 

0.40 (1.02) 

0 (0) 

47 (119) 

26.32 (66.85) 

27.0 (68.6) 

109.07 

(703.68) 

0.95 (2.41) 

0 (0) 

46 (117) 

26.88 (68.28) 

28,0 (71.1) 

99.16 

(639.74) 

0.91 (2.31) 

0(0) 

62 (157) 

31.38 (79.71) 

34.5 (87.6) 

293. 1 1 

(1891.03) 

1.56 (3.96) 

0 (0) 

49 (124) 

36.28 (92.15) 

39.0 (99.1) 

94.89 

(612.19) 

0.9Q (2.29) 


Andropogon virginicus^ 
(Broomsedge) 

Stalk 

Leaf 

Melinis minutifJora 
(Molasses grass) 
Stalk 
Leaf 

Pennisetum setaceum 
(Fountaingrass) 

Stalk 

Leaf 
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Table 2 Minimum and maximum, mean, median, variance, and standard error of the mean for stalk and leaf 
surface area-to~volume ratios of three major grass fuels in Hawaii 


Species 

Min 

Ma.x 

Mean 

Median 

Variance 

Standard 
error of 
the mean 

A ndropogon virginicus 
(Broomsedge) 

Stalk 

533 

(17) 

lift (Ijcm) 

6000 (197) 1064 (35) 

906 (30) 

541,661 (583) 

74 (2.4) 

Leal 

4000 

(131) 

12000 (394) 

7036 (231) 

6000 (197) 

3,575,256 (3848) 

189 (6.2) 

Mel in is n i inut iflo ra 
(Molasses grass) 

Stalk 

471 

(15) 

1200 (39) 

676 (22) 

658 (22) 

16,714 (18) 

13 (0.4) 

Leaf 

3000 

(98) 

8000 (262) 

4536 (149) 

4400 (144) 

972,947 (1047) 

99 (3.2) 

Pennisetum setaceum 
(Fountaingrass) 

Stalk 

326 

(11) 

962 (32) 

590 (19) 

583 (19) 

15,677 (17) 

13 (0.4) 

Leaf 

2182 

(72) 

6000 (197) 

3447 (113) 

3429 (112) 

534,373 (575) 

73 (2.4) 


Table ^ — Bed depth and surface area-to-volume 
ratios, by species - 


Table - Weightirig factors ^ leaf and stalk, and for live and 
dead categories, by species 


Species 

Weighting factors 

l.eaf 

Stalk 

Live 

Dead 

Andropogon virginicus 

0.6634 

0.3366 

0.4463 

0.5537 

Melinis minut {flora 

0.3465 

0.6535 

0.3252 

0.6748 

Pennise turn setaceum 

0.4168 

0.5832 

0.3405 

0.6595 


Species 

Bed depth 

Surface 

area-to- 

volume 

ratio 


Feet (m) 

///> (Hem) 

Andropogon 

virginicus 

1.624 (0.495) 

5026 (165) 

Melinis 
minut (flora 

2.210 (0.674) 

2013 (66) 

Pennisetum 

setaceum 

2.785 (0.849) 

1781 (58) 


Table 5 — Minimum and maximum, mean, median, variance, and standard error of the mean fuel weight on base subplots, and the corresponding 
loading, by species 


Min 1 





Standard error 



Ma.K 

Mean 

Median 

Variance 

of the mean 

Loading 


(kg) Lbffd (kginf) 


A ndropogon virginicus 
Mel inis m in utiflo ra 
Pennisetum setaceum 


0.300 (0.136) 
1.138 (0.516) 
0.514 (0.233) 


0.723 (0.328) 
4.159 (1.887) 
2.344 (1.063) 


0.427 (0.194) 
2.332 (1.058) 
1.376 (0.624) 


0.404 (0.183) 
2.119 (0.961) 
1.312 (0.595) 


0.0169 (0.0035) 
0.9987 (0.2055) 
0.3198 (0.0658) 


0.041 (0.019) 0.124 (0.605) 

0.316 (0.143) 0.677 (3.305) 

0.179 (0.081) 0.400 (1.953) 


univariate distribution. The means 
correspond to Equation 3. The me- 
dians also are estimated with respect 
to the {c/xj. 

The refinements to surface area-to- 
volume that we estimated required 
more effort than was originally 
planned. In view of the nonlinear 
dependence of Rothermel’s rate of 
spread model on however, the 
additional work is justified. 


RESULTS 

We found pronounced differences 
in the fuel characteristics of leaves and 
stalks, particularly for surface area- 
to-volume ratios {tables 1 and 2). We 
weighted the mean values for leaf and 
stalk according to the ovendry weights 
of grab samples of the plant that were 
segregated into leaf and stalk masses, 
and also into live and dead masses 


{table 3). Estimates of bed depth and 
surface area-to-volume ratio in tables 
1 and 2, weighted with the leaf and 
stalk proportions in table i, produce 
the weighted means in table 4. On the 
basis of our study, they represent the 
estimates of bed depth and surface 
area-to-volume ratio that are appro- 
priate for a fuel model of each. 

Finally, we estimated fuel loading 
from the base subplot samples {table 







Table 6 — Mmimiim and maximum, mean, median, variance, and standard error of the mean fuel weight for all subplots, and the corresponding 
loading, by species 


Species 

Min 

Max 

Mean 

Median 

Variance 

Standard error 
of the mean 

Loading 

Andropogon virginicus 

0.143 (0.065) 

Lb (kg) 

0.524 (0.238) 0.273 (0.124) 

0.253 (0.115) 

0.01 11 (0.0023) 

0.033 (0.015) 

Lb /ft" (kg/m-) 

0.079 (0.386) 

Melinis minutiflora 

0.458 (0.208) 

3.481 (1.579) 

1.721 (0.781) 

1.786 (0.810) 

0.8807 (0.1812) 

0.297 (0.135) 

0.500 (2.441) 

Pennisetum setaceum 

0.347 (0.157) 

0.966 (0.438) 

0.616 (0.279) 

0.600 (0.272) 

0.0487 (0.0100) 

0.070 (0.032) 

0.179 (0.874) 


Table 1—Dead and live weights and corre- 
sponding loadings, by species 


Species 

Weight 

Loading 

Andropogon 

virginicus 

LJy (kg) 

Ljhift- (kg/m-) 

Dead 

0.151 (0.068) 

0.044 (0.215) 

Live 

Melinis 

minutiflora 

0.122 (0.055) 

0.035 (0.171) 

Dead 

1.161 (0.527) 

0.337 (1.645) 

Live 

Pennisetum 

setaceum 

0.560 (0.254) 

0.163 (0.796) 

Dead 

0.406 (0.184) 

0.118 (0.576) 

Live 

0.210 (0.095) 

0.061 (0.298) 


5) . The sample size for each species is 
10, which represents the number of 

base subplots per transect. The mean 
values in the table cannot be used 
without recognizing that the samples 
on which they are based (base sub- 
plots) had relative maximum volumes. 
To compensate for this bias, we use 
estimates of fuel volume the “nonbase” 
subplots had relative to the base sub- 
plots. Let Wj be the loading estimate 
and W’i^j the airdry weight of the base 
subplot material at stake i on the 
transect; let {Vjj} be the volume pro- 
portions, relative to the base subplot, 
on the nonbase subplots at stake i. 
Then 

3 

Wi = (1 + Z Vij)WiJ4 

j=i 

The {wj}, averaged in the usual man- 
ner, yield the loading statistics {table 

6) . We tried two approaches to 
apportioning the loadings among live 
and dead categories. In the first, we 
used observations of percent green, 
and in the second, we used the 


proportions-by-weight data for live 
and dead (table i). The analysis of 
simultaneous observed estimates of 
percent green and proportions-by- 
weight indicated that the observed 
estimates are consistently lower than 
the latter, by an average of 8.9 
percent. Differences between the two 
approaches were analyzed by using 
the Wilcoxon signed rank test.'^ For 
eight differences, we observed a 
Wilcoxon signed rank statistic, w = 2. 
The corresponding 2-tailed signifi- 
cance probability was 0.0234, which 
we consider significant. Ocular esti- 
mates of live and dead loading pro- 
portions are apparently not easy to 
make, without appreciable statistical 
error. (By error, we mean the differ- 
ence between the proportions of live 
vegetation estimated (a) by observa- 
tion, and (b) by weight measurements.) 

The dead and live loadings reported 
were calculated using the proportions- 
by-weight (table 7). 


CONCLUSIONS 

Characteristics of the fuels we 
studied differ considerably from those 
of the standard fuel models. The 
surface area-to-volume ratio measured 
in broomsedge is 2.51 times as large as 
that in the National Fire-Danger 
Rating System model for perennial 
grass, and 1.44 times as large as the 
largest area-to-volume ratio of the 
models described by Albini.^5 [jge of 
the reciprocal harmonic mean of 
thickness measurements to estimate a, 
which we justify, gives larger values 
than the reciprocal arithmetic mean. 
The dead fuel loading we measured in 
molasses grass is larger than any of 


the other loadings for dead 1-hour 
fuels in the standard models. Calcula- 
tions by Burgan,’^ using, Hilo, Hawaii 
weather data, indicate that the fuel 
characteristics reported in this study 
yield significantly higher fire danger 
indices than those produced by Model 
L of the 1978 NFDRS.i? Further 
study is needed to determine the effect 
of the new estimates on criteria for the 
adjective fire danger classes. This will 
require a comparison of the empirical 
distributions of the respective fuel 
models. 
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